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SUMMARY 


An aerodynamic investigation of a parabolic body of revolution was 
conducted at a Mach number of 1.92 with and without an annular supersonic 
jet exhausting from the base. Measurements with the jet inoperative were 
made of lift, drag, pitching moment, radial and longitudinal pressure 
distributions, and base pressures. With the jet in operation, measure- 
ments were made of the pressures over the rear of the body with the pri- 
mary variables being angle of attack, ratio of jet velocity to free- 
stream velocity, and ratio of jet pressure to stream pressure. 

The results with the jet inoperative showed that the radial pres- 
sures over the body varied appreciably from the distribution generally 
employed in most approximate theories. The linearized solutions for lift, 
pitching moment, and center of pressure gave relatively poor predictions 
of the experimental results. An analysis of several theoretical methods 
for calculating pressure distribution and wave drag showed that some 
methods gave results in considerable disagreement with experimental values. 

Maximum effects of the jet were obtained at the lower ratio of jet 
velocity to stream velocity and the highest ratio of jet pressure to 
stream pressure. These effects amounted to a slight decrease in fore- 
drag, a reduction in lift, and a shift of center of pressure in a 
destabilizing direction. 


INTRODUCTION 


Aerodynamic investigations at supersonic speeds of bodies of 
revolution simulating those containing jet— propulsion units have almost 
entirely neglected the effects of the jet flow upon the flow over the 
rear of the body. An experimental subsonic investigation of the effects 


-'■Supersedes recently declassified NACA Research Memorandum L9K09 
by Eugene S. Love, 1950* 
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of the jet upon the aerodynamic characteristics of the aggregate 
A-5 missile (body plus four equally spaced tail surfaces) was conducted 
in Germany in 19^0. (See ref. 1.) The results of these tests showed 
the jet to cause (l) an increase of as much as 100 percent in the normal 
forces at small angles of attack, (2) a shift of center of pressure to 
the rear by an average of about 0.5 maximum body diameter, and (3) 321 
increase of drag of approximately 70 percent. Other investigations, both 
subsonic and supersonic, of jet effects upon the flow over bodies were 
conducted with the A-b missile at an angle of attack of 0°. (See ref. 2.) 
The results of subsonic drag tests were in general agreement with those 
found in tests of the A-5 missile. The jet caused an increase of drag of 
as much as 80 percent. The results of the supersonic tests showed a maxi- 
mum decrease of drag of 18 percent. 

In most instances, the best aerodynamic design of bodies housing jet 
units entails a certain degree of boattailing; that is, convergence of 
the body surface as it approaches the jet exit. Boattailing such that 
the diameters of the jet exit and of the exterior body surface become 
equal would probably favor greater jet effects upon the flow over the 
rear of the body than would other geometric conditions; therefore, it was 
chosen as the geometric condition to be employed in the present investi- 
gation. 

The primary purpose of the investigation was to determine the effects 
of an annular supersonic jet exhausting from the base of a parabolic body 
of revolution upon the flow over the rear of the body. It was necessary 
to obtain first the aerodynamic characteristics of the body without the 
jet. Therefore, comprehensive force and pressure-distribution measure- 
ments were made of the basic jet model body. Similar but not as exhaustive 
auxiliary tests were conducted on a parabolic body (same body family but 
larger thickness ratio) initially employed during bench tests of small 
annular nozzles developed for use in the present investigation. All tests 
were conducted in the Langley 9- inch supersonic tunnel at a Mach number 
of 1.92. The Reynolds numbers for the tests were 2.5I x 10 b for the jet 
model body and 2.A7 X 10 6 for the auxiliary body. 


SYMBOLS 


A, B 


a 


constants in equation of parabola defining body shapes (eq. ( 1 ) ) 
angle of attack 


P = \ 
C D 


- 1 


, /Drag \ 

total drag coefficient -5 

M^tiax/ 
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C Dmin minimum drag coefficient 


C D b 

CDp 

^Df 

C DW 


base drag coefficient (Pb 


Smaxy 


minimum foredrag coefficient ^CD min - (^ c 'Db) a -QoJ 

j Sy \ 

skin- friction drag coefficient Cfs „ 

V llam Smax / 


wave-drag coefficient 


Cf lam 


skin-friction coefficient for laminar flow on a flat plate 
1-328 

Vr 


C L 


total lift coefficient 


/ Lift 

\qSmax 


Cl 


weighted unit lift 


Cjj, pitching-moment coefficient Momen^about reference point 

1 «c 


dmax maximum body diameter 

6 lift density 

e apex half-angle of body 

l body length 

l c cutoff body length 

1^ completed body length (tip to tip) 

M free-stream Mach number 

M(b es design jet Mach number based on area ratio 


k 
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n 

9 

Pa 

Pe 

Ps 

P°m 

P 

Pb 

*1 

<1 

P 

r 

R 

s b 

Sm 

Smax 

^w 

t 


angle of local- surface inclination with respect to axis of 
symmetry- 

radial angle 

atmospheric pressure 

static pressure of jet at jet exit 
stream pressure or pressure of ambient air 
pressure in model stilling chamber 
pressure increment 

pressure coefficient 

base pressure coefficient 
lifting-pressure coefficient 

dynamic pressure 

density of fluid 
radius of body 

Reynolds number referred to l c 
base area 

mean cross-sectional area for body of length l c 
maximum frontal area 
wetted area of body of length 

thickness ratio 
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v^ c volume for cutoff body length 

Vj^ volume for completed body length 

V undisturbed stream velocity 

Vj velocity of jet 

x longitudinal coordinate 

Subscripts : 

max maximum 

I>=0 value at zero lift 

APPARATUS AND TESTS 
Wind Tunnel and Model Installation 


The Langley 9- inch supersonic tunnel is a closed-return direct-drive 
tunnel in which the pressure and humidity of the enclosed air may be con- 
trolled. Throughout the tests the quantity of water vapor in the tunnel 
air was kept at sufficiently low values so that negligible effects on 
the flow from condensation were present in the supersonic nozzle. The 
test Mach number is varied by means of interchangeable nozzle blocks 
forming test sections approximately 9 inches square. A schlieren optical 
system provides qualitative visual flow observations. Eleven fine-mesh 
turbulence-damping screens are installed' in the settling chamber ahead 
of the nozzles. 

Figure 1 shows the general installation for tests of the jet model. 
Pressure within the model stilling chamber was varied by means of manually 
controlled valves installed ahead of the juncture of the incoming- air- 
supply line with the flexible-air-supply line. Force and pressure- 
distribution measurements of the models with jet inoperative employed 
the same model support system with the air-supply system removed. The 
scales used are self-balancing beam scales and measure three components, 
in a horizontal plane, of the total forces on the model and support 
system. 
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Description of Models 

All models were constructed of mild steel, were highly polished, 
and, except for a special pressure-distribution model, were mounted on 
slender, hollow sting supports which, for the jet model, served also 
as an air- supply conduit. The surface contours of the models were deter- 
mined by revolving about its chord a parabolic arc obtained from the gen- 
eral parabolic equation 


r 


Ax - Bx 2 


( 1 ) 


In this equation the constants A and B can be easily obtained for 
desired values of maximum diameter, base (or jet-exit) diameter, and 
thickness ratio. (See appendix.) 


Three separate models were constructed with a surface contour given 

hy 


r = 0.l827x - O.OI85 Ax 2 (2) 

The designations assigned these models were: model 1-J, the basic jet 

model with two interchangeable tail sections containing jet nozzles 
of Maes = 2.11 (nozzle 1) and M^ es = 3*19 (nozzle 2); model 1-F, the 
model employed in the force tests; and model 1— P, a special pressure- 
distribution model constructed in two halves about a meridian plane and 
containing 63 pressure orifices located in one-half of the model along 
three meridians, 0°, 45°, and 90°, with 21 orifices similarly spaced 
along each meridian. 

The auxiliary model tested had a surface contour given by 

r = 0.2460 k - 0.02647* 2 (3) 

This model was designated model 2. 

Except for model 1-P, the bases of all the models were hollow or 
open, as for the case of a jet exit. Special plugs were made to fill 
the annular base openings of models 1— F and 2 flush with the body ends 
for use in tests of these models with a simulated solid or closed base. 
Photographs of models 2, 1-J, and 1-P are shown in figure 2. Model 1-F 
has been excluded since its external appearance is no different from 
that of model 1-J. 
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The following table gives the pertinent geometric parameters of 
the models: 


Parameter 

'Model 1 

Model 2 


7.719 

7.607 

h> in 

9.854 

9.293 

t 

-0.09135 

0.1230 

e, deg . 

IO.36 

13.83 

Vj , cu in. . 

3.113 

4.857 

Vj , cu in 

3.3^6 

5.080 

Sy, sq in 

16.339 

20.330 

Snp sq in. . 

0.4036 

O.6385 

St>, sq in 

O.2923 

0.3526 

W aq in 

O.6365 

1.0272 

^max-’ in 

0.9002 

1.1436 


Development of Annular Nozzles 

Numerous- bench tests were conducted to determine suitable shapes 
and sizes of annular nozzles that might be constructed in the tail sec- 
tion of model 1-J. Design of a theoretically shock- free annular super- 
sonic nozzle contour of such small size was not attempted in view of the 
analytical complications, boundary- layer effects, and the difficulty of 
machining to the desired accuracy a curving, internal contour of such 
small radii. Nozzle 2 (Mdes = 3*19) represented the best attempt at con- 
struction without prohibitive surface imperfections of a nozzle with a 
curving contour to give the jet a flow direction at the exit similar to 
that of nozzle 1 (Mdes = 2.11). In spite of extreme care, small imper- 
fections in the surface contour of this nozzle could be detected. Because 
of insufficient pressure of the air-supply facility, conclusive bench 
tests of nozzle 2 could not be made. The higher ratios of p e /p s (ratio 

of jet static pressure to pressure of ambient air) obtainable for nozzle 1 
allowed reasonably conclusive bench tests of this nozzle. Diametrical 
surveys at the nozzle exits were conducted by means of a 0.010- inch total- 
pressure tube mounted in a micrometer traversing arrangement. The total- 
pressure tube measured pressures on a plane perpendicular to the nozzle 
center line and just beyond the nozzle lip. Static pressure within the 
jet was measured by means of an orifice vented to the nozzle just inside 
the lip. The Mach number distribution across the nozzle exits was cal- 
culated from these pressures with the assumption of negligible effects 
due to the slight difference in longitudinal positions of the static- 
and total-pressure measurements and the assumption that the static pres- 
sure across the jet was constant. For the values of p e /p s of the bench 
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tests, a conically shaped nozzle was found to give the most uniform dis- 
tribution at the jet exit for a design Mach number of 2.11. Figures 3(a) 
and 3(b) show the results of nozzle surveys from bench tests and from 
tests conducted in a similar manner with the use of the tunnel as a partly 
evacuated container for the model to obtain larger values of p e /p s * The 
surveys using the tunnel as a vacuum chamber (p s ~ 0.6p a in fig. 3(a) 
and p s ~ 0.5p a in fig. 3(b)) show a marked improvement in the distri- 
bution for nozzle 2 and a slight lessening of the "hump" in the distri- 
bution curve for nozzle 1. The marked improvement in the distribution 
for nozzle 2 is apparently a result of the decrease in the pressure rise 
across the shock originating at the lip of the nozzle and reflected by 
the sting surface, and a decrease in the boundary- layer buildup caused 
by back pressure which in turn tends to eliminate compressions in the 
flow within the nozzle. In the bench tests of certain of the annular 
nozzles of M^ es = 3 or greater (p s = p a ), the large pressure rise across 
the shock from the lip caused a thickening of the boundary layer near the 
lip of the outer nozzle surface and a region of reverse flow that extended 
a considerable distance away from the inner (sting) surface. For the 
cases for which reverse flow could not be detected, the results indicated 
that the large adverse pressure rise across the shock caused a rapid thick- 
ening of the boundary layer along the sting surface ahead of the point of 
reflection of the shock. The adverse pressure gradients and the thick- 
ening boundary layers probably caused compressions in the flow ahead of 
the shock and a resulting rapid drop in velocity at the outer and, par- 
ticularly, the inner diametrical stations. (See fig. 3(b).) With the 
decrease in external pressure (p s « 0 . 5p a ) the pressure rise across the 
shock from the lip decreases; therefore, the adverse effects would also 
be expected to decrease in a manner similar to that indicated in fig- 
ure 3(b) . By similar reasoning, a satisfactory Mach number distribution 
would be expected at- the exit of nozzle 2 in the tunnel tests at M = 1.92 
for which p s » 0.lAp a . 

Pressure measurements indicated that the best position for the 
orifice measuring the pressure in the model stilling chamber was 

that shown in figure 1. Thermocouple measurements showed that the 
temperature of the air in the model stilling chamber varied very little 

from storage— tank air temperature. Values of the reference pressure p n 

m 

for the jet tests were measured by means of a large Bourdon gage. An 
open- tube manometer, used in conjunction with this gage, served as a 
constant check of the pressure gage and supplied values of p Qm for 

pressures less than atmospheric. Figure k shows the calibration curves 
for each nozzle with the tunnel in operation. Although the values 
Tojq were intended to serve only as accurate reference pressures, 

figures A and 5 show that they have some quantitative value as well. 

The values of Mach number calculated from values of p e / p and presented 
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in figure 4 for nozzles 1 and 2 compare favorably with the average values 
of the Mach number distributions of figures 3( a ) and. 3(b) .> respectively. 
In addition, figure 5 shows that the thrust of nozzle 1 obtained at two 


values of 


■P°m/ -Pa 


by calculations based upon p and the Mach number 

°m 


distribution checks closely the thrust measured by strain-gage apparatus 
during the bench tests. 


Tests 


All tests were conducted through an angle-of-attack range of approxi- 
mately ±5°. Mirrors approximately l/l6 inch square were flush mounted in 
the bodies near the base as a part of the optical angle-of-attack system. 
Force tests and base-pressure measurements of models 1-F and 2 were made 
with base open and base closed for three longitudinal positions of the 
models. These were body base even with, l/2 inch ahead of, and 1 inch 
ahead of the end of the sting windshield. All drag values were corrected 
for the buoyancy effect due to the difference between free-stream pressure 
and the pressure within the box enclosing the sting windshield and balance. 
Radial and longitudinal surface-pressure measurements were made with 


model 1-P at meridian intervals of 22± along every meridian from 0° 
o ^ 

to 112 ^ and at l8o° (0° to l80° represents angle-of-attack plane). With 


the jet in operation, the base of the model was 1 inch from the end of 
the sting windshield. The primary variables of the jet tests were a, 
Pe/Ps> anc ^ Vj/V. For the measurements of the jet effects upon the pres- 


sures over the rear of the body, the tubes were installed as shown in 
the inset in figure 1. Previous investigations showed that the lead 
tubes in such an arrangement had no measurable effect upon the pressures 
over the body along a meridian l8o° opposite. All schlieren photographs 
were taken with the knife edge horizontal. 


Precision of Data 

The estimated probable errors in the aerodynamic quantities are 
included in the following table. The value of ±0.08° given for angle 
of attack is a result of error in the initial referencing of the model 
bodies with respect to stream direction. The value of +0.01° is the 
error that might be incurred in relative angle-of-attack readings for 
a given test. The values for Cl, Cl, and C m apply only to the 

results obtained from the mechanical scales. 
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C L 

C D 


M 

CL, 

deg 

R 

P 

Initial 

Relative 

±0.0004 

±0.0004 

±0.0018 

±0.01 

±0.08 

±0.01 

±20,000 

±0.002 


Comparison of the actual ordinates of the model bodies with the values 
obtained from equations (2) and (3) showed the body dimensions to be 
accurate, with one exception, within ±0.002 inch. This exception, the 
tail section of model 1-J containing nozzle 1, had gradually increasing 
small errors in the radii of the body from a point approximately 0.3 inch 
from the base rearward. The maximum error in radius (at the body base) 
amounted to ±0.008 inch. The effects of this lesser degree of boat- 
tailing will be shown in the results. The meridian planes and rotational 
angles for the radial pressure distributions were accurate within ±2°. 

The Bourdon gage for measuring p^ gave readings accurate within 

±0.2 pound per square inch. 


RESULTS AND DISCUSSION 
Jet Inoperative 


Force tests .— Figure 6 shows the aerodynamic characteristics of 
model 1-F for the three longitudinal positions of the model in relation 
to the forward tip of the sting windshield. Corresponding schlieren 
photographs for these and two additional positions are shown in figure 7. 
In like manner, the aerodynamic characteristics of model 2 and schlieren 
photographs at two longitudinal positions are shown in figures 8 and 9 , 
respectively. Except for the zero longitudinal position, all schlieren 
photographs in figures 7 and 9 were taken at zero angle of attack. Values 
of Cm in figures 6 and 8 are for moments taken about the point of 
maximum diameter. 

The results of base— pressure measurements with varying longitudinal 
position and angle of attack indicated much the same effects from the 
presence of the sting support and windshield and from angle pf attack 
as discussed in references 3, 4, and 5; specifically, the base pressures 
vary appreciably with angle of attack, and the body undergoing test must 
be mounted on a long, slender sting support if base pressures simulating 
free-flight values are to be obtained. The results of the force tests 
are given in table I. 


NACA TN 3709 


11 


The condition of the base of the bodies, open or closed, had little 
or no consistent effect upon the results except a slight increase in the 
base drag for the closed condition. The values of base drag at the zero 
longitudinal position of the models are to serve only in establishing 
the magnitude of the foredrag and not as accurate measurements of the 
base drag since the proximity of the sting windshield to the body base 
would affect the base pressures. The effects on the lift and moment 
curves from the flow impinging upon the exposed sting at the higher 
angles of attack are shown in figures 6 and 8. These effects increase 
as the exposed sting area increases and cause marked nonlinearities in 
the curves. The longitudinal position of the models apparently had 
little effect on the minimum foredrag coefficient although close 

examination of the schlieren photographs in figures 7 and 9 shows a 
lessening of the laminar separation near the base of the body with 
increasing distance between the body base and the tip of the sting 
windshield. 

Values of drag coefficient due to skin friction Cp^ (laminar flow 

was observed over the entire body) were calculated for the test Reynolds 
number. These values and their approximate percentage of the foredrag 
are presented in the following table: 


Model 

c Df 

Percent of Cjjp, 

1— F 

0.0216 

30 

2 

.0167 

14 


References 6, 7, and 8 have pointed out independently that at least 
to first order the limiting value for the lift- curve slope of very slender 
bodies of revolution at small a is (expressed in radians and based 
upon Sb) 





and that the center— of— pressure location in relation to the nose of the 
body is 


Center of pressure = 1 — [ — u 

W 


(5) 
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From equations (4) and (5) it follows that the slope of the pitching- 
moment curve, with moments taken about the nose of the body, is 



(6) 


The values calculated from equations (4), (5)> and (6), expressed in 
degrees and referred to Smax.’ are presented in the following table and 

compared with the experimental values (in parentheses) obtained at the 
1-inch longitudinal position. The experimental values of lift- curve 
and moment-curve slopes given in table I include support interference 
effects and aerodynamic tares on the exposed sting. However, the experi- 
mental slope values are for zero lift, and pressure measurements along 
a 1-inch length of the exposed sting from the body base have shown the 
lifting forces upon the sting to be negligible within an angle -of- attack 
range of +2°. Furthermore, the effects of the presence of the sting and 
windshield upon the body lifting forces would be expected to be least 
at the 1-inch position. 


Model 

10 ^ 

da 

C 0 p 0 

(&iam. from nose) 

about nose 

da 


0.0160 

3.26 

- 0.00610 

1-F 

(.0290) 

( 1 . 28 ) 

(-. 00431 ) 


.0120 

5.49 

-.00975 

2 

(.0338) 

( 1 . 39 ) 

(-. 00705 ) 


All the theoretical values are relatively poor predictions of the experi- 
mental results. Part of the failure of equations (4), (5), and (6) to 
predict values in reasonable agreement with experimental values is prob- 
ably caused by the use of the geometrical value of base area. Refer- 
ence 3 has shown that the calculated pressures over the rear of a body 
of revolution with boattailing, as given by the method of characteristics, 
are in excellent agreement with experimental pressures if the pressure 
calculations are performed along the streamline of separated boundary 
layer. This would seem to indicate that the geometrical value of base 
area in equations (4), (5), and (6) should be replaced by an area deter- 
mined by the diameter between the separated streamlines at the body base. 
Measurements of this "diameter of separation were made from enlarged 
schlieren photographs of models 1-F and 2, each at the 1-inch longitudinal 
position and with open base. An area of 0-347 square inch was obtained 
for model 1-F and an area of 0.474 square inch for model 2. The values 
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calculated from equations (4), (5); and- (6) and referred to these areas 
are presented in the following table and compared with the experimental 
values (in parentheses) : 


Model 


c.p. 

(diam. from nose) 

d-Cm , , 

about nose 

da 

da 

1 — TT 

0.0190 

1.40 

-0.00310 

-i-j? 

(.0290) 

(1.28) 

(-.00431) 

p 

.0161 

2.31 

-.00559 


(.0338) 

(1.39) 

(-.00705) 


Although these values are an improvement upon the previous theoretical 
values, they are still rather poor predictions. 

Pressure distributions .- The results of radial-pressure-distribution 
measurements are presented in figure 10 for model 1-P and in fig- 
nr p ii for model 2. Longitudinal pressure distributions are presented 
in figures 12 and 13 for model 1-P and in figure 14 for model 2. Radial 
pressures for model 2 are given in figure 15 . Although the results for 
model 2 are secondary to those for model 1-P, they tend to indicate that 
certain phenomena observed in the pressure distributions of both bodies 
apparently hold for slender pointed bodies of revolution in general. 

First, figures 10, 11, 13(b), 14(c), and 15 show that the pressures along 
the 90° meridian at a = 0° do not remain relatively unchanged with 
angle of attack, a simplifying assumption often employed in approximate 
theories for computing the aerodynamic characteristics of conical bodies 
and pointed bodies of revolution. In fact, at certain horizontal stations 
the pressure at the 90° meridian varies as much as or more than at any 
other meridian. Second, the radial pressure distribution at any longi- 
tudinal station varied appreciably from the usually assumed cosine dis- 
tribution, especially for the low-pressure half of the body at longi- 
tudinal stations ahead of the maximum thickness and for the entire 
circumference at stations behind the maximum thickness. The radial pres- 
sures at stations behind the maximum thickness behave in much the same 
manrifer as observed in tests of slender cylinders in yaw. . (See ref. 9*) 
Third, for longitudinal stations ahead of the maximum thickness there 
appears to be a radial point for each longitudinal station at which the 
pressure remains almost constant with angle of attack. (See fig. 15 > 
in particular.) The locus of these radial points does not follow a 
meridian but tends to move away from the 90° meridian in the direction 
of the high-pressure half of the body as the distance from the nose 
increases. For model 1-P, this shift was from 0 « 80° at station 0-0°° 
to 0 » 48° at station O.606. For model 2, the shift was from 0 » 7° 
at station 0.283 to 0 « 6l° at station 0-573. The fourth phenomenon 
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observed, was "the consistent appearance of the hump, not predicted, by 
potential theory, in the longitudinal- pressure-distribution curves. For 
model 1-P, this characteristic occurred near the 0.3 station, and for 
the a = 0° condition amounted to a noticeable discontinuity in the 
curve. For model 2, it occurred in the vicinity of the 0-5 to 0.6 sta- 
tions. This phenomenon has been present in the results of other tests 
of slender pointed bodies of revolution (refs. 5 and 10, for example) 
and, except for cones, is apparently characteristic of slender pointed 
bodies of revolution in general. Although cross velocities in the vicin- 
ity of the 90° meridian would not be expected to affect the lift, their, 
inclusion would, nevertheless, be expected to reduce the pressures at 
the 90° meridian, possibly of the magnitude observed in the experimental 
results. Also, when the experimental longitudinal pressure gradients 
in the vicinicy of the 90° meridian are found to be relatively large, 
the experimental tangential pressure gradients are found to be of the 
same order of magnitude. A tangential gradient of such magnitude would 
be expected to have important bearing upon separation effects. 

Figure l6 illustrates the method by which the pressure coefficient 
at any point on the body is converted to the lifting -pres sure 
coefficient P^. The equation, including second— order terms, is 

Pj = P(cos 0 cos T) cos a + sin T) sin a) (7) 

All values of P for model 1-P were converted to P 7 by means of equa- 

tion (7) . With Smax as the reference area, the total lift coefficient 
would be given by (see diagram in upper half of fig. 17) 

Ct = — — — ^ f P 7 r d0 — - — dx (8) 

L ^max Jo Jo 1 cos n 


Plots of cos q against 0 were made for each horizontal station 

and graphically integrated from 0 to 2 jt. This gave what might be 
termed the lift density 5 at each station. Therefore, the weighted 
unit lift for each station x is 
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From the linear first-approximation theory of reference 7 an l the first- 
order approximation theory of reference 11, a solution for the lift dis- 
tribution over the body from equation (4-) would apparently give an 
acceptable first-order prediction. In the dimensions of equation (9) 
and for small values of a, this solution may be expressed as 


(=l)x - ‘“' r *(I) x (10) 


In figure 17 the values of (c^) x from equations (9) and. (10) are plotted 

against horizontal station, in inches, for values of a of 2. 50° and 5 • 00® . 
This gives a graphical representation of the lift distribution over the 
body. Values of Cl and C ra determined by integration of the experi- 
mental curves and the theoretical values (in parentheses) are presented 
in the following table : 


a, 

deg 

°L 

CL 

m 

about nose 

c.p. 

(diam. from nose) 

2.50 

0.0717 

(.0400) 

-0.00733 

(-.01525) 

0.88 

(3.26) 

5.00 

O -3 
0 ro 

-.0230 

(.03050) 

1.34 

(3.26) 


Since the tests of model 1-P most nearly duplicate in tunnel posi- 
tion the tests of model 1-F at the 1-inch longitudinal sting station, 
a reasonable check of the pressure-distribution results should be realized 
in a comparison with the force-data results for the 1-inch sting station 
(open base). In an effort to eliminate as much as possible the effects 
of the sting support, values based upon the slopes of the curves at zero 
lift were determined from table I. These values are included in the fol- 
lowing table. Of course, the use of constant slopes will mask the effects 
of separation at the rear of the body . 


CL, 

deg 

9 L 

Cm 

about nose 

c.p. 

(diam. from nose) 

Source 

2.50 

0.0725 

-0.0108 

1.28 

Constant slope 

5.00 

.1450 

-.0216 

1.28 

Constant slope 

2.50 

.0700 

-.00766 

.94 

Average— curve value 

5.00 

.1450 

-.0146 

.86 

Average— curve value 


16 


NACA TN 3709 


Model 1-F was also tested as shown in figure 18. The results are given 
in figure 19 and in the following table: 


a, 

deg 

°L 

about nose 

c.p. 

(diam. from nose) 

2.50 

0.0860 

-0.0246 

2.45 

5.00 

.1831 

-.0532 

2.49 


These values do not show close agreement with the values obtained from 
model 1-P, but it is interesting to note that they agree closely with 
the results of the tests of model 1-F, open base, at the zero longi- 
tudinal position; thus, an appreciable effect from the flush sting- 
windshield arrangement is indicated. As might be expected, the only 
close check of the values obtained through integration of the pressure 
distributions lies in a comparison at a = 2.50 0 with the average-curve 
values from tests of model 1-F at the 1-inch longitudinal sting station. 
On this basis, the results of the pressure-distribution integrations may 
be considered reliable. The shift of center of pressure with angle of 
attack, as determined from the pressure distributions, would correspond 
to effects that might be expected from separation of the flow at the 
rear of the body. 

The results of the pressure-distribution investigations give some 
insight into the causes of the higher experimental values of dCp/da, 
as compared with theoretical values, noted in these and other tests of 
slender pointed bodies of revolution. Figure 17 shows that, for body 
stations behind the station at which the theoretical local lift would 
be a maximum, the experimental values of local lift exhibit a marked 
increase over the theoretical values. The station of initial gain in 
the experimental local lift also appears to be in the region where the 
experimental radial pressures begin to deviate appreciably from a cosine 
distribution. It is seen that the observed noncosine distributions of 
radial pressures act in such a way as to give more lift over the body 
than would cosine distributions. Figure 17 also shows the expected 
reduction of antilift forces by separation in the region over the rear 
of the body where re compress ion would otherwise take place. 

Analysis of theoretical methods for prediction of longitudinal 
pressu re distribution .— The equations for several methods for predicting 
the pressure distribution over slender pointed bodies of revolution were 
calculated in terms of equation (1). Reference 8 gives a solution 
termed the "rigorous linearized first-order solutipn" for an extremely 
sharp nose body of revolution for M > \f2. This yields 


P » — 2(A2 - 6ABx + 6B 2 x 2 )log(A - Bx) 


(ID 
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which was obtained, in reference 8 from a previously derived equation for 
the pressure coefficient for compressible flow, given herein as 


P = -A 2 + l6ABx - 22B 2 x 2 + 2(A2 - 6ABx + 6B 2 x 2 ) |log 2 - log(0A - 0Bx)J + 


0(A3 - 3A2Bx + 3AB 2 x2 _ b3x3) 


( 12 ) 


By the method of reference 12, but with axes at the nose of the body, 
the general equation was found to be 


P = 2 


.|(6AB — 9B 2 X ) \j x 2 — p 2 r 2 + [a 2 — 6ABx + 3B^(2x 2 + B 2 r2)J cosh 


(13) 


Reference 13 gives a solution for P that is identical with equation (13) 
except for one additional term. This solution is 

P = Value from equation (13) - (A 2 - 4ABx + 4B 2 x 2 ) (14) 

Reference 10 has given a solution based upon the small-disturbance theory 
and requiring a step-by-step numerical integration. For simplicity this 
method will be expressed herein merely as 


p -zIE< r ’Mj] 

i=l 


(15) 


where i and n represent the initial and nth integration station, 
respectively. A 25— point calculation was made. The final method 
employed was the characteristics method of references 14 and 15 with the 
assumption of potential flow. 
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Equations (ll) to (15) were applied to the body shapes of models 1 
and 2. The characteristics method was applied to model 1 only. The 
results of these pressure— distribution calculations are presented in 
figure 20. The corresponding wave drag coefficients from inte- 

gration of the pressure curves are given in the following table: 


Method 

Model 1 

Model 2 

Equation (ll) 
Equation (12) 
Equation (13) 
Equation (l4) 
Equation (15) 
Characteristics 
method 

0.0253 

.0784 

.0822 

.0740 

.0746 

.0624 

0.0356 

.1498 

.1548 

.1302 

.1309 


Integration of a curve determined by the experimental points of 
figure 20 gives for model 1 a pressure-drag coefficient of 0.0542 and 
for model 2 a drag of 0.11. (This value is very approximate because of 
lack of sufficient points.) It should be noted that no correction has 
been applied to the experimental points since surveys have shown only 
minor variation in the static-pressure distribution in the wind-tunnel 
test section. Therefore, any correction applied to the pressure drag 
would be minor. The effects of separation upon the experimental pressure- 
drag coefficient would not oppose the indication that all these theories 
predict too large a wave drag. As might be expected, the characteristics 
method shows best agreement with experiment. Although it gives only a 
fair prediction, equation (l4), the Lighthill solution, is the best of 
the approximate theories and gives a slightly better prediction of the 
experimental results than does the laborious small-disturbance method of 
equation (15). Indications are that equation (ll) should not be applied. 

If the values of the experimental pressure drag are assumed to 
approach the wave drag, the addition of the laminar skin— friction drag 
should give a value that checks closely the measured foredrag. The 
comparison is given in the following table. The corresponding values 
of the foredrags are from table I, 1-inch position, open base. 


Model 

Experimental pressure 

Foredrag from 

drag plus Cpj. 

force teats 

1 

O.O758 

0.0743 

2 

.1267 

.1104 
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Jet in Operation 

Figure 21 presents schlieren photographs of the jet model with 
and without tubes to the pressure orifices installed as shown in the 
upper left-hand corner of figure 1. As previously stated, pressure 
measurements with no jet throughout the angle— of-attack range showed 
no effect upon the body pressures from the presence of or disturbances 
produced by the pressure lead tubes on the side of the body opposite 
the pressure orifices. The surveys and calibrations of the two jet 
nozzles indicated reasonable values of the Mach numbers for the two 
jet nozzles to be approximately 2.10 and 3.O5. For a free-stream Mach 
number of 1.92, these values would represent ratios of jet velocity to 
stream velocity Vj/V of approximately 1.05 and 1.24, respectively. 

Figure 22 shows the pressure change at each orifice location due 
to jet effects with varying jet pressure and angle of attack. Also 
included are the hysteresis effects (at the a = 0° and a = 2.50°, 
upper surface, for both velocity ratios) dependent upon whether the 
particular test was made with increasing or decreasing jet pressures. 

For both velocity ratios the major effects upon the pressures over the 
rear of the body occurred at a = 0° and a = 2.50°, upper surface, 
and were confined to the orifices nearest the jet exit. At these a con- 
ditions and a velocity ratio of 1.05, the body pressures increased posi- 
tively as the jet pressure approached and exceeded stream pressure, the 
greatest change occurring immediately after the over-pressure condition 

Pe 

“ ^ "Weis reached. At the sane a, conditions and a velocity ratio 

-uS 

of 1* 24, the body pressures showed a very slight decrease at the under— 

Pe 

pressure conditions — <1 and a slight increase at the over-pressure 

Ps 

conditions. At a = 5*00° the jet had no significant effect at either 
velocity ratio. Figure 23 gives the percentage change in body pressures 
due to the jet at a = 0° and a = 2.50°, upper surface, for several 
values of p e /p s . At the top of the figure the differences in the basic 
pressure distributions over the rear of model 1-P and the two nozzle tail 
sections (no jet) of model 1-J are presented for a = 0°. These differences 
in pressures appear to be the effects of previously mentioned very small 
measured differences in body ordinates. The zero-percent datum lines of 
the plots of jet effects represent the pressures with no jet. The major 
effects of the jet upon the body pressures are confined to approximately 
5 percent of the body length (from the base) for a velocity ratio of 1.24 
and to approximately 10 percent of the body length for a velocity ratio 
of I.05. For similar pressure ratios p e /p s the effect of the jet upon 
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the body pressures in these regions is much greater for the lower velocity 
ratio. At a = 2 . 5 O 0 , upper surface, there is a positive increase in jet 
effects over the a = 0° condition. This might be expected in view of 
the greater separation of the flow from the upper surface at angle of 
attack that would favor pressure effects from the jet to be felt farther 
forward along the body and to a greater degree. 


Typical schlieren photographs of the jet in operation at a = 0° 
are shown in figure 2k for the two velocity ratios and, whenever possible, 
for similar pressure ratios. Photographs at the lower velocity ratio 
were taken with the lead tubes installed; therefore, for comparison with 
the photographs at the higher velocity ratio which were taken without 
the tubes installed, the simple pattern of disturbances present in the 
jet-inoperative photograph should be ignored. 


As the supersonic jet begins to flow, there is a decrease in thick- 
ness of the laminar boundary layer at the rear of the body. This is 
particularly true at the higher velocity ratio. Up to a value of p e /p s 

equal to approximately 0 - 5 , the shock pattern within the jet and at the 
nozzle lip is much the same for both velocity ratios. A very prominent 
lambda shock is noted at the jet outer boundary immediately rearward of 
the nozzle lip. The shock pattern within the jet follows closely the 
expected phenomena for under-pressure jets calling for the presence of 
a shock to offset the higher pressure outside the jet boundary . From 


0.5 to 1 . 0 , the lambda shock at the jet outer boundary tends toward 
Ps 

a plain shock whose apparent point of origin at the jet boundary lies 
slightly downstream of the nozzle lip. The shock pattern within the jet 
continues along the pattern for under-pressure jets except that two shocks 
are now observed within the jet of higher velocity ratio. It is possible 
that this may he due to slight imperfections in the nozzle surface; but, 
in view of the fact that a similar phenomenon, though not as strong, may 
be observed at slightly higher pressure ratios at the lower velocity, 
ratio, it appears more probable that this phenomenon is characteristic 


of the flow from annular nozzles of this type. 


At a value of — 

P c 


1 , 


the iet boundaries are parallel to the axis of symmetry of the jet. 

P 

Parallel flow would not be expected to occur exactly at = 1 because 

Ps 

the stream pressure would be somewhat different from the pressure behind 
the shock emanating from the vicinity of the nozzle lip. Also, the pres- 
ence of the shocks within the parallel jet at ~ 1 may be attributed 

Ps 

to the necessary change in flow inclination from a direction that is 
away from the axis of symmetry Immediately ahead of the exit at the 
nozzle lip to a direction parallel to the axis of symmetry beyond the 
exit. At the higher values of p e /ps separation occurs at the rear of 
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the body; this separation is more pronounced at the lower velocity ratio. 

An expansion of the jet flow occuirs at the nozzle lip for the over-pressure 

p 

conditions — > 1 and is followed by shocks within the jet. (See 
Ps 

fig. 25-) With increasing over pressures of the jet, the shock from the 
outer boundary of the jet near the jet exit reverts more and more to a 
lambda shock at the higher velocity ratio, whereas the corresponding 
shock at the lower velocity ratio continues as a clearly defined plain 
shock. Indications from the pressure measurements and the schlieren 
photographs are that the induced velocities imparted to the semidead 
air in the region where the flow is separated from the rear of the body 
are greater for the case of the higher velocity ratio. Apparently, these 
induction effects tend to counteract the back-pressure effects over the 
body from the related internal (jet) and external flow characteristics. 


The effects of the jet upon the foredrag of the body (a = 0°) were 
determined from the measured body pressures. The results are tabulated 
in the following table as percentage change in : 


Vj/’ 

Wpvx 

1.05 

1.24 

0.4 

0.34 

— 

1.0 

-.21 

0.72 

1.5 

-1.01 

0 

1.8 

- 1.62 

-.43 

3-5 

-3.31 

— — 


The maximum effect of the jet upon the lift and pitching moment occurred 
at a = 2 . 50 ° and, as in the case for the foredrag, at the lower veloc- 
ity ratio, 1 . 05 ? and highest pressure ratio, 3 » 5 * An approximation of 
this effect is entered on the curve representing lift distribution in 
figure 17 . The results of such a change in lift distribution would be 
a 13 . 7 -percent decrease in the overall lift of the body and a destabilizing 
shift of center of pressure from 0.88 diameter rearward of the nose to a 
point 0.74 diameter ahead of the nose. Although the foregoing determina- 
tion of jet -effects has assumed the interference effects of the sting and 
windshield on the pressure measurements to be the same for the jet-on 
and jet-off conditions, the differences in the effects are expected to 
be small. Thus, the foregoing values should give an insight into the 
magnitude of some of the effects that might be expected from an exhausting 
annular jet. 
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CONCLUSIONS 


An aerodynamic investigation at a Mach number of 1.92 of a parabolic 
body of revolution with and without a supersonic annular jet exhausting 
from the base indicates the following conclusions for the case with the 
jet inoperative: 

1. The condition of the base of the body, hollow or closed, has 
little and no consistent effect upon the aerodynamic characteristics of 
the body. 

2. The simplified linearized solutions for lift-curve slope, 
pitching-moment- curve slope, and center-of -pressure location give rela- 
tively poor predictions of the experimental results. 

3. Experimental radial pressure distributions show a marked devia- 
tion of lifting pressures from the theoretical cosine distribution, and, 
contrary to the simplifying assumption of most approximate theories, the 
pressures on the sides of the body ( 90 ° from angle -of- attack plane) vary 
appreciably with angle of attack. These discrepancies appear to be the 
result of separation effects and of the failure of the theories to include 
effects of cross velocities which may be important. 

4. The method of characteristics for axial symmetry gives a reasonable 
overall prediction of the actual pressure distribution over the body. How- 
ever, a "hump," not predicted by potential theory, is found in the experi- 
mental longitudinal-pressure-distribution curve at forward body stations. 
This phenomenon appears to be characteristic of slender pointed bodies of 
revolution in general. 

5 . The Lighthill solution appears to be the best of the linearized 
solutions investigated for prediction of pressure distribution over slender 
bodies of revolution. 

For the case with the jet in operation, the following conclusions 
are indicated: 

6. Pressures over the rear of the body show little effect from the 
jet until the jet pressure ratio approaches and exceeds the value for 
parallel flow of the jet. 

7 . The effects of the jet upon the body pressures are reduced as 
the ratio of jet velocity to free-stream velocity increases. 

8 . The greatest effects of the jet upon the lift and pitching moment 
occur at 2 . 5 O 0 angle of attack and almost completely disappear as the 
angle of attack is increased to 5 * 00 °* 


NACA TN 3709 


23 


9. Maximum effects of- the jet are obtained at the lower ratio of 
jet velocity to stream velocity, 1.05, and the highest ratio of jet pres- 
sure to stream pressure, 3*5* These effects amount to a 3* 3-percent 
reduction in foredrag at 0° angle of attack and, at 2.5O 0 angle of attack, 
a 13. 7-percent reduction in lift and a destabilizing shift of center of 
pressure from 0.88 body diameter rearward of the nose to a point approxi- 
mately 0.7^ diameter ahead of the node . 

10. Indications are that, for higher ratios of jet velocity to stream 
velocity than achieved in the present investigation, the induction effects 
of the jet upon the flow over the rear of a body of this type would become 
as important as the back-pressure effects. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 25, 19^9. 
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APPENDIX 
Body Parameters 


The general equation for the shape of the bodies is given as 


r = Ax - Bx 2 


(Al) 


The constants A and B may be determined simply as follows. At maxi- 
mum thickness and in nondimensional form 


r max _ A 2 
Z kBl 


(A2) 


Also 


A _ base 

B = 2— (A3) 

l 

If r , rwooo, and the thickness ratio are assigned fixed values, 
max' oas s / a 

the constants A and B are readily obtained by combining equations (A 2) 
and (A3) . Examination of the basic equation (Al) shows that the con- 
stant A is dimensionless and is equal to twice the thickness ratio. 
However, the constant B is not dimensionless. Therefore, any calcula- 
tions employing equation (Al) with the numerical values replacing the 
constants A and B, such as computations of pressure distributions or 
drag, must be carried through with the same dimension units (feet or 
inches) used in calculating the values of A and B. 

The following equations apply to the family of bodies determined 
from equation (Al) : 


Volume 



(Ak) 


A 2 1 2 _ abz3 + B 2 Z^\ 
V 3 ' 2 5 j 


(A5) 


Wetted surface area = it 




(A6) 
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TABLE I.- SUMMARY OF FORCE-TEST RESULTS FOR 
MODEIB 1-F AMD 2 WITH JET INOPERATIVE 


Model 

Distance from 
"body base to 
f orward end of 
sting windshield 
(in.) 

Base 

condition 

(— ) 

\d<x y L=0 

^ 4=o 

atout Vx 

C ^mln 

M a=o° 

c.p. 

(diam. from nose) 


0 

Open 

0.0334 

0.0121 

0.0899 

0.0142 

2.37 


0 

Closed 

.0330 

.0124 

.0879 

.0156 

2.25 

1-F 

.5 

Open 

.0324 

.0126 

.0940 

.0214 

2.14 


• 5 

Closed 

.0323 

.0102 

.0955 

.0233 

2.77 


1.0 

Open 

.0290 

.0142 

.0960 

.0217 

1.28 


1.0 

Closed 

.0292 

.0140 

.0960 

.0235 

1.45 


0 

Open 

.0392 

.0112 

.1370 

.0217 

2.16 


0 

Closed 

.0394 

.0110 

.1373 

.0220 

2.21 

2 

.5 

Open 

.0427 

.0105 

.1400 

.0241 

2.43 


• 5 

Closed 

.0429 

.0105 

.1400 

.0252 

2.43 


1.0 

Open 

.0338 

.OI 36 

.1410 

.0306 

1.39 


1.0 

1 

Closed 

.0343 

.0130 

.1407 

.0309 

1.94 
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(a) Models 2 and 1— J. 



(b) Side view of model 1— P. 



(c) 


Auxiliary view of tube exits, model 1— P. 


Figure 2.- Models. 



L-630^1.1 
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Taper from throat -3.12° 



Figure 3«- Results of surveys of nozzle exits. 
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Ps 



Figure 3«- Concluded. 
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Figure A.- Calibration of jet nozzles in tunnel. 
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Figure 5.- Thrust of nozzle 1. 
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(a) Base to windshield, is 0 inch. . 

Figure 6.- Aerodynamic characteristics of model 1-F . 
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(b) Base to windshield is 0.5 inch. 
Figure 6.- Continued. 
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(c) Base to windshield is 1.0 inch. 
Figure 6.- Concluded. 
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(a) 0 inch. 



(b) 0.5 inch. 



(c) 1 inch. 



L-6j0k2 


Figure 7«- Schlieren photographs shoving the effect upon the flow over 
the rear of model 1-F with varying distance between body base and tip 
of sting windshield. 
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(a) Base to windshield is 0 inch. 

Figure 8.- Aerodynamic characteristics of model 2. 


NACA TN 3709 


39 


<3 

*\ 

,NI 


V 


r 










ju\ 


c 




A 


L 





T // 



a 



p 


P 


/V — ** 


r 

1 

fl 









1 








Base { 
open ' 

r oC/_ 

□C/77 about d/r/Qx 

A Cr\ 




if 











N. t/ 

F/ags , base c/osed 

— ^ — 











A 


X 













r 



) 














• 














r 




P 

X 











A 

r' 


Jp 










[ 

X 



// 

// 

u 













-Jf 














// 

// 































































\NACAP^ 

1 L 


OC deg 


(b) Base to windshield is 0.5 inch. 
Figure 8.- Continued. 
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(c) Base to windshield is 1.0 inch. 
Figure 8.- Concluded. 
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(a) 0 inch. 



(b) 0.5 inch. 



L- 630li.3 


Figure 9 .- Schlieren photographs showing the effect upon the flow over 
the rear of model 2 with varying distance between body base and tip of 
sting windshield. 
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Figure 10.- Radial pressure distributions for model 1-P. 
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Figure 10.- Concluded. 
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(a) Body stations ahead of maximum thickness station. 
Figure 11.- Radial pressure distributions for model 2. 
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(b) Body stations behind maximum thickness station. 


Figure 11.- Concluded. 
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Figure 12.- Variation of longitudinal pressure distribution with angle 
of attack at two meridians of model 1-P. 
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Figure 13 .- Longitudinal pressure distribution of model 1-P at a = 0° 
and variation of 90 ° meridian pressures with angle of attack. 
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Figure Ik.- Longitudinal pressure distributions for model 2. 
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Figure 15 •- Radial pressure distributions for model 2. (Examples of 
radial point of zero pressure variation with angle of attack.) 
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Figure l6.- Conversion of normal pressures to lifting pressures. 
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Figure 17- - Graphical representation of the lift distribution over 

model 1-P. 
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Figure l8.- Model 1-F with flush sting windshield. 
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Figure 19 •- Aerodynamic characteristics of model 1-F; surface of sting 

windshield flush with body base. 
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(b)Mode/ 2. 


Figure 20.- Comparison of theoretical and experimental pressure distribu- 
tions for a = 0°. (Numbers in squares indicate equation numbers in 
text . ) 
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Figure 21.- Schlieren photographs of phenomena at jet exit without and 

V i 

with pressure lead tubes installed; — = 1.24. 
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Figure 22.- Effects of jet upon pressures at rear body stations with varying 

angle of attack and jet pressures. 
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Figure 22.- Concluded 
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Figure 23.- Effects of jet upon rear body pressures. 
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(a) No jet. 




(a) No jet. 



Figure 2k. - Schlieren photographs of jet in operation at two ratios of 
jet velocity to free-stream velocity and with varying ratios of jet 
pressure to free-stream pressure, a = 0°. 
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Figure 24. - Continued 
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Figure 24.- Concluded. 
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(a) Schlieren photograph of flow phenomena. 



(h) Sketch of flow phenomena. 
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Figure 25.— Phenomena in region of jet exit; = 1.05* — = 3«1°* 
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